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Abstract 

Phi meson production in pion-baryon and baryon-baryon interactions near 
the production threshold is studied in a boson-exchange model. Parameters 
in this model are either taken from the Bonn model for the nucleon-nucleon 
potential or determined from the measured phi meson decay width and the 
experimental data for the reactions n~p^(f)n and pp^ppcj). The isospin- 
averaged cross sections for the reactions TrN—xf)N, TrA^<f)N, NN^NN(j), 
which are needed in transport models to study phi meson production in heavy- 
ion collisions at subthreshold energies, have been evaluated. For the reactions 
NA^NN(p and AA^NN(f>, which become singular when the exchanged pion 
is on-shell, the Peierls method is used to calculate their cross sections. We 
find that the cross section for the reaction nN^KjjN is much larger than the 
other cross sections. 
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The production of phi mesons, which are pure ss states, in hadronic interactions is usually 
suppressed due to the Okubo-Zweig-Iizuka (OZI) rule ]1J. An enhanced production of phi 
mesons in ultra-relativistic heavy-ion collisions has thus been suggested as a possible signal 
for the deconfmement phase transition of hadrons to a quark-gluon plasma 0, where the 
strangeness production rate has been shown to be order-of-magnitude larger than in a free 
hadronic gas. However, the phi meson properties might be changed in hot dense matter due 
to the partial restoration of chiral symmetry. Indeed, studies based on QCD sum rules |||J , 
hadronic models including vacuum polarization effects [f|, and the hidden gauge theory || 
have all shown that the phi meson mass is reduced in medium. The change of phi meson 
properties can be studied in dilepton spectra 0] and may also be inferred from the phi meson 
yield in heavy- ion collisions ||. 

Phi meson production from heavy-ion collisions has been studied in CERN experiments 
from the dimuon invariant mass spectra |J. It has been found that the ratio for <p/u> in S+U 
collisions at 200 GeV/nucleon is enhanced by about a factor of 3 compared to that from the 
proton-proton and p-U interactions at the same energy. The enhancement can probably be 
accounted for by assuming that a quark-gluon plasma has been formed in the collisions. On 
the other hand, it can also be explained by hadronic scenarios P,r0|,n|. In particular, it has 



been pointed out in Ref. |§ that the reactions KA — > (pN and KK — > <pp become important 
if one takes into account the medium effects on hadron masses in hot dense matter. 

Phi mesons have also been measured recently in heavy-ion collisions at AGS energies 
from the invariant mass distribution of kaon-antikaon pairs [[12[], and their mass and width 



are found to be consistent with that in free space. This is not surprising as these kaon- 
antikaon pairs are from the decay of phi mesons at freeze out when their properties are the 
same as in free space. If a phi meson decays in medium, the resulting kaon and antikaon 
would interact with nucleons, so their invariant mass is modified and can no longer be used 
to reconstruct the phi meson. 

Phi meson production from heavy-ion collisions at SIS/GSI energies are being studied 
by the FOPI collaboration also through the K + K~ invariant mass distribution | 13fl - Their 



preliminary results indicate that the <fi/K~ ratio is about 10%, and is thus similar to that 
observed at AGS energies. This is quite surprising as the incident energies at GSI and AGS 
differ by order of magnitude. In a previous study by us |n| , the phi meson yield in central 
Ni+Ni collisions at 2 GeV/nucleon has been studied via its decay into dileptons. Including 
phi meson production only from the kaon-antikaon annihilation channel, the (f)/K~ ratio 
was found to be about 1%, which is much below the value of 10% from the preliminary data 
of the FOPI collaboration Jl3|]. It is thus necessary to also consider phi meson production 



from other processes, such as the baryon-baryon and pion-baryon interactions. 

There are a number of advantages to work at heavy-ion collisions at SIS energies, which 
are below the threshold for strange particle production from nucleon-nucleon interactions 
in free space. First, it has been shown that subthreshold particle production in heavy ion 
collisions is sensitive to medium modifications of hadron properties |15|-p^. Second, at 
these energies the reaction dynamics is relatively simple as the colliding system consists 
mainly of nucleons, delta resonances, and pions. In addition to kaon-antikaon annihilation 
as considered in Ref. |TJ]], phi mesons can be produced from the pion-nucleon, pion-delta, 



nucleon-nucleon, nucleon-delta, and delta-delta collisions. The cross sections for these reac- 
tions are thus needed as inputs in the transport model. Unfortunately, these cross sections 
are not known empirically near the threshold, which are relevant for phi meson produc- 
tion in heavy ion collisions at subthreshold energies. In this letter, we shall introduce a 
boson-exchange model to calculate these cross sections. Similar models have been used pre- 
viously to study pion |L8|], kaon fl9|-plj|, e ta meson p2fl , and dilepton production in 



nucleon-nucleon interactions in the same energy region. 

Specifically, we shall calculate the following cross sections: irN-^N<f) , 7iA^N(j) , 
NN^-NNcj) , NA^NN<p , and AA—>-NN<j) . In principle, one needs for heavy ion col- 
lisions the inclusive cross sections such as NN — > (f)X instead of the exclusive cross sections 
such as NN^NNcj) . Since the energies at SIS are below the phi meson production threshold 
in the nucleon-nucleon interaction, we expect that exclusive reactions such as NN^NNcf) , 
which have lower thresholds, are more important than the ones with pions also in the final 



state, which require additional energies and are thus suppressed. We note that the boson- 
exchange model can be generalized to the case with one or two pions in the final state using 
the approximation of resonance dominance, namely, replacing one or both nucleons in the 
final state by the delta resonance. Such studies will be pursued in the future. 

The Feynman diagrams for nN^-Nip and 7rA— >N(p are shown in Fig. 1, while those 
for NN^NN<p ,NA^NN<p , and AA^NN<f) are shown in Fig. 2, Fig. 3, and Fig. 4, 
respectively. For baryon-baryon interactions there is also the exchange diagram as the final 
nucleons are identical. The Lagrangians for nNN, pNN, pNA f24|, and 7rp0 |25| interactions 
are well-known, and they are given by 

C-kNN = -^^^7 5 7 M ^ • ^vf, (1) 
C-pNN = -9pNNi>l tl r^/j ■ p p - ^^-^a^Tij; ' ( d nPv ~ d »Pit)> ( 2 ) 

c pNA = i^^^fr ■ (d^ v - d v p,) + h.c, (3) 

Ulp 

= ^e^d^dafo ■ ir. (4) 

In the above, ip is the nucleon field with mass m^; ip^ is the Rarita-Schwinger field for the 
spin-3/2 A resonance; and n, p p , M are meson fields with masses m n , m p , m^, respectively. 
The isospin matrices, Dirac 7-matrices and Levi-Civita tensor are denoted by r, 7 M and 
£ fivap ^ reS p ec tively. 

The matrix element for the reaction TrN^N(f) includes both the vector and tensor con- 
tributions and is given by 

WpNNj-Kp4> 



M 



fpNNfirptp 



u(p f )Yu{pi)\ e af 3 lu q a e (3 k' y 

u(p f )a^k p u(p i )]e a ( hv q a e p k' ) '. (5) 



2m N m ( f ) (t — m 2 

In the above, u and u are Dirac spinors of the initial and final nucleons with momenta Pi and 
Pf, respectively. The momentum and polarization of the phi meson are denoted by q and e, 
respectively, while t is the square of the four-momentum k of the exchanged rho meson. 
The differential cross section is then 



do 1 1 , , rl9 ,„ s 

Tt = 6to^ |M| • (6) 

where and \fs are, respectively, the pion three-momentum and the total energy in the 
center-of-mass frame of the pion-nucleon system. In the above equation, the sum over final 
spins and average over initial spins are implied. We note that in evaluating \M\ 2 in the above 
equation and also that for other reactions discussed below, the program Form (Version 1.0) 
|26| has been extensively used. 

The coupling constants (f n NN, QpNN, f P NN, I p na) and cut-off parameters in C^nn, 
£ p nn, and £ p na are taken from the Bonn one-boson-exchange model (Model II) as listed in 
Table B.l of Ref. [24]. From the measured width T^ np pa 0.6 MeV, the coupling constant 



fnp4> ~ 1-04 is determined. In addition to the form factor at the pNN vertex which is taken 
from the Bonn model, we also introduce a monopole form factor with a cut-off parameter 
A^ p ^ at the npcf) vertex. By fitting to the available experimental data for the reaction 
7i~p — > <pn [pjj, we obtain A£ ^ = 1.2 GeV. The comparison of the calculated cross section 
(solid curve) with the data (open circles) is shown in Fig. 5. We see that the dominant 
contribution comes from the tensor pNN coupling (long-dashed curve) as that from the 
vector pNN coupling (short-dashed curve) and the interference term (dot-dashed curve) 
are small. In the same figure, we also show the parameterization introduced by Sibirtsev 
|28| (dotted curve). Although both are similar at energies where data exist, they differ 
appreciably near the threshold. 

Having determined the cut-off parameter, we can calculate the isospin- averaged cross 
sections for 7tN—>N(j) and irA^Ncf) without introducing further free parameters. In our 
study the A particle is treated by the usual Rarita-Schwinger formalism, with its spin 
projection operator given by 

The results are shown in Fig. 6. We see that the cross section for TcN—^Ncj) (solid curve) 
is about a factor of five larger than that of 7rA— >N(p (dotted curve), mainly because of the 
tensor pNN coupling in the former reaction. 



The matrix element for NN^NN(p includes both the direct and exchange terms and is 
given by 



M = M d + M e , 



(8) 



where 



M d = 



firNN9pNNfirp(j> 



m 



2d 



ifizNNfpNNfizpif) 



2m N m 7r m (t> (q 2 d - ml)(q 2 d - m 2 ) 



u{pz)lz>iidu{pi) u{pi)a p pu{p 2 



tpvafiq^^vqid 



T pi/af3 



and 



qp^v(q2d) a q2d, 



fnNN9pNNfnp^> 



(9) 



mvm<t>{qle ~ ™l){q\e ~ ™ 2 P ) 

ifnNNfpNNfnp<t> 

2m N m- K m 4> (q 2 e - m\){q\ e - m 2 p ) 
e^qM^Ul- 



M(p 4 )75^1eM(Pl)J [U(p3h f3 u(p 2 ) 

u{pi)l<oiieu{pi) u(p 3 )a p/3 u(p 2 ) 



tpvapq^^vq^e 



(10) 



In the above, p\ and p 2 are the four-momenta of two initial nucleons; p 3 and p 4 are those 
of final nucleons; and q and e are, respectively the momentum and polarization of the phi 
meson. The momenta of the exchanged pion and rho meson are denoted, respectively, by 
qid — P3 - Pi and q 2d = p 4 - p 2 for the direct term and by q le — p 4 - p 1 and q 2e = p 3 - p 2 
for the exchange term. 

For the reaction NN — > NN<p, in addition to the usual form factors at the 7rNN and 
pNN vertices taken from the Bonn potential model, two additional monopole from factors 
are introduced at the npcfi vertex as both pion and rho meson are off-shell. The one associated 
with the rho meson is taken to be the same as in the reaction tcN — > N<p. The other one 
associated with the pion introduces another cut-off parameter A£ ^, and it can be determined 
by fitting to the cross section for the reaction pp — > ppcfi. Unfortunately, no experimental 
data are available for this reaction near the threshold where the boson-exchange model is 
expected to be appropriate. Nevertheless, we choose A^ p<}!) =0.95 GeV to fit the experimental 
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Piab ~ 10 GeV p9y30|, which is the lowest beam energy with experimental data 
available. The resulting cross section is shown in Fig. 7 by the solid curve together with 
the experimental data (open circle). It should be noted that A£ ^ and are both in the 
order of 1 GeV, similar to typical values for cut-off parameters in the Bonn potential model. 

The cross section for pp — > ppep has also been studied by Sibirtsev based on an one- 
pion-exchange model. In this study, both the off-shell feature of the pion and the interference 
between the direct and exchange diagrams are neglected, and the cross section for pp — > ppcf) 
can thus be expressed in terms of the n~p — > <pn cross section, which is taken from the 
empirically measured one. This approach, first introduced in Ref. ||31|| , has been used earlier 



for calculating kaon production cross section in the nucleon-nucleon interaction p9[pi| . The 
results of Ref. |28j (dotted curve) are seen to be significantly larger than our results, which 
includes both the off-shell and interference effects. Experiments on phi meson production 
from the proton-proton collision near the threshold are being carried out at SATURN |32| . 
These data will be very useful in checking the validity of our model. 

The extension of the model to phi production in nucleon-delta (A r A) and delta-delta 
(AA) interactions is straightforward, except for the complication arising from the fact that 
at the 7riVA vertex the energy momentum conservation allows the pion to go on-shell. A 
pole thus develops in the pion propagator at some region of the phase space, which then 



leads to a singular cross section. Similar singularities appear also in iVA — * NAK [21] and 



NA — > NNt] p2[ . In both Ref. and Ref. [22], a complex pion self-energy has been 



introduced to remove the singularity. Since the pion self-energy is density dependent, the 
resulting cross section becomes density dependent, and diverges when the density approaches 
zero. This treatment also introduces certain model dependence as different pion self-energies 



have been used. In Ref. p2"|, the pion self-energy is related to the pion-nucleon scattering 



cross section, while in Ref. pTJ it is calculated in a simple A-hole model. 



In this work, this singularity is removed by using the so-called Peierls method [|33 
which has also been used very recently in the treatment of the singularity in the reaction 
— > evW + associated with the physics of muon collider [34}]. The basic idea of the 



Peierls method is to introduce a complex four-momentum for the resonance in the initial 
state to account for its finite lifetime. As a result, the four-momentum of the exchanged 
particle acquires an imaginary part through the energy momentum conservation. In our 
case, the pion propagator in NA—>-NN<j) and AA—>-NN<j) is modified as follows: 

1 1 

(Pa - Pn) 2 ~ ml (p A - p N ) 2 -ml- im A {EA ~f N) T A ' 

where p^ and pa are the four momenta of the final nucleon and the initial delta resonance, 
respectively; En and Ea are their energies in the center-of-mass frame; and tua is the 
delta mass. The singularity at {jpA — Pn) 2 — m\ is seen being removed by the finite delta 
width Y a ~ 120 MeV. The resulting cross sections are shown in Fig. 8. Near threshold 
the NA—*NN<j) cross section (dashed curve) is about a factor of five larger than that of 
NN^>NN(j) (solid curve), while the A A— >NN<fi cross section (dotted curve) is only slightly 
larger than the latter. 

In summary, we have calculated the phi meson production cross sections in pion-baryon 
and baryon-baryon interactions based on a boson-exchange model. Most parameters in the 
model are taken from the Bonn model for the nucleon-nucleon potential. Two additional cut- 
off parameters are introduced at the irp<j) vertex, and are determined by fitting to available 
experimental data. We have calculated the isospin averaged cross sections for ttN-^-NcJ) , 
7iA^N(p , NN->NN<j) , NA^NN<p , and AA^NN<p without introducing any further 
parameters. It is found that the cross section for irN—>-N<f) is much larger than the other 
cross cross sections for phi meson production. These cross sections will be useful for the 
transport model study of phi meson production in heavy-ion collisions at SIS energies, and 
such a study will be reported elsewhere f35fl . Also, to study possible medium modifications 
of the phi meson properties, in particular the reduction of its mass, the measurement of 
dilepton spectra will be very useful fli"i| , |36|| . Experimentally this will be carried out by the 
HADES collaboration for heavy- ion collisions at SIS energies ||37|| . A detailed transport 
model study of dilepton production including phi meson production from pion-baryon and 
baryon-baryon interactions and from other processes is also in progress. 
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FIGURES 





FIG. 1. Feynman diagrams for TrN^Ncp and irA^Ncj) 




FIG. 2. Feynman diagrams for NN^NN<ft 





FIG. 3. Feynman diagrams for NA^NNfi . 




FIG. 4. Feynman diagrams for AA^NN(j) 
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FIG. 5. The cross section for it p — > (ftn 
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FIG. 6. Isospin-averaged cross sections for ivN—>N<t> and 7rA— >iV</> . 
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FIG. 8. Isospin-averaged cross sections for NN^NNcf) , NA^NNcf) and AA^NNcj) . 



